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We have calculated the potentials of the heavy (charmed or bottomed) pseudoscalar mesons up to 
0(e 2 ) with the heavy meson chiral perturbation theory. We take into account the contributions from 
the football, triangle, box, and crossed diagrams with the 2(f) exchange and one-loop corrections to 
the contact terms. We notice that the total 2<^>-exchange potential alone is attractive in the small 
momentum region in the channel BB 1 ^ 1 , _B s _B s /_0 , or _B_B S /_1//2 , while repulsive in the channel 
BB I ~° . Hopefully the analytical chiral structures of the potentials may be useful in the extrapolation 
of the heavy meson interaction from lattice QCD simulation. 
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I. INTRODUCTION 

Since the discovery of X(3872) pQ, many charmonium-like and bottomonium-like states such as JT(3940) [2], 
X(4160) [3], et al. have been observed in the past decade. The charmonium-like state X(3872) was first observed 
by the Belle Collaboration in the exclusive decay process B ± — > K ± 7r + 7r~ J/ip. Last year the Belle Collaboration 
observed two charged bottomonium-like resonances Z(,(10610) and Z;,(10650) in the hidden-bottom decay channels 
7r ± T(nS') (n = 1, 2, 3) and n ± h b (mP) (to = 1, 2) of T(55) 0]. 

Some of these new states including the above two charged Zb states do not fit into the conventional quark model 
framework. Various theoretical approaches including the lattice QCD the QCD sum rule [5] and the quark 
model [7] have been employed to interpret the underlying structure of these new states. Despite huge experimental 
and theoretical efforts, the nature of some of these exotic states is still elusive. 

For example, the interpretation of X(3872) remains challenging since the discovery in 2003. One popular speculation 
is that X(3872) is a molecular state composed of a pair of heavy mesons [STUOj. Similarly, the two charged Zf,(10610) 
and Zb( 10650) states are proposed as the BB* and B* B* molecule states within the one boson exchange (OBE) 
framework [Til H2] • 

Besides the charmonium-like and bottomonium-like states, the possible existence of some molecular candidates 
composed of BB mesons and DD mesons is also very interesting. If the attractive interaction is strong enough 
between the heavy meson pair, this kind of states may exist. Their behavior will be very similar to the deuteron 
which is composed of two nucleons. There have been some investigation of these interesting states within the OBE 
model. 

However, the interaction potential derived from the OBE model contains several phenomenological coupling con- 
stants and cutoff parameters, which should in principle be extracted through fitting to the experimental data. Unfortu- 
nately, there is not much experimental information on the strong interaction between the light meson and heavy meson. 
It will be very desirable to derive the strong interaction between the heavy meson pair with a model-independent 
approach. Especially many new states such as X(3872) and the two Zb states lie very close to the threshold. Within 
these very loosely bound systems, the long-range pion exchange force should play an important role. Therefore the 
chiral perturbation theory provides a natural framework to investigate the heavy meson strong interaction. In this 
work, we shall derive the heavy pseudoscalar meson potential order by order. 

Chiral perturbation theory (%PT) is a model-independent tool to study the chiral dynamics of heavy hadrons. 
Heavy hadron xPT is frequently used for the system made up with a single heavy hadron and light pseudoscalar 
mesons because of its explicit power counting |13H17| . The scattering matrix can be expanded order by order in the 
small parameter e = p/A x , where p represents either the momentum of the light pseudoscalar mesons or the residual 
momentum of the heavy hadrons in the nonrelativistic limit, while A x represents either the scale of chiral symmetry 
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breaking or the mass of heavy hadrons. The power counting guarantees that one can just calculate some limited 
Feynman diagrams and obtain the scattering matrix at the certain order. 

Weinberg developed a new formalism and first extended the chiral perturbation theory to the two nucleon system 
[T51 IT!?] . Since his pioneering work, the modern nuclear force has been built upon the chiral effective field theory [2"01 - 
130] . Such a formalism was used extensively to investigate the various few-nucleon observables such as the partial- wave 
analysis, few-nucleon scattering, reaction etc. 

As pointed out by Weinberg [TB] US] , the power counting of the two nucleon scattering matrix is broken by the 
double poles of the heavy hadrons in some 2-particle-reducible (2PR) Feynman diagrams. Let's illustrate this point 
with the box diagram in Fig. [I] The Feynman amplitude can be written as 



a L io + P o + l£ _ l o + P o +l£ x ~ l J aL io + P o + l£ _ l0 + P o + l£ J ai > W 

where we omit the parts relevant to the pion which preserves the power counting. We will focus the integral with I , 
and work it out by closing the 1° contour integral in the lower half-plane 



1 ' : ' ' l° + P° + ie -1° + P° + ie ~ P° + ie ~ pz/(2m N )+ie' 



The naive power counting predicts that X should be 0(1/\P\). But it is 0(rriN /\P\ 2 ) from Eq. (|2j)! X is actually 
enhanced by a large factor ttin/\P\ compared to the naive power counting prediction. 



P + l 
P 
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P 



FIG. 1: The box diagram. The solid line represents the nucleon, and the dashed line represents the pion. 

With Weinberg's formalism, we do not directly calculate the scattering matrix of the few hadrons perturbatively 
with the heavy hadron %PT due to the 2PR diagrams. Instead, we focus on the potential. In the derivation of 
the hadron-hadron potential, one takes into account the 2-particle-irreducible (2PI) parts of the Feynman diagrams 
only and calculate the potential of few hadrons perturbatively with the correct power counting. Afterwards, one 
can obtain the scattering matrix with the potential by solving the nonperturbative equations such as Schrodingcr 
equations, Lippmann-Schwinger equations, and so on. The 2PR contributions will be recovered when solving the 
nonperturbative equations. 

The reliable hadron-hadron potential is a necessary input for getting the scattering amplitude or phase shift of 
the hadrons. It is also essential to explore the existence of the heavy hadron molecules. For example, the binding 
energy or size of the molecular states can be obtained from the potentials of the hadrons by solving the Schrodingcr 
or Lippmann-Schwinger equations. 

In this work, we shall use Weinberg's formalism to derive the BB potentials in four independent channels up to 
I-loop level with heavy meson chiral perturbation theory (HMxPT). We include the heavy vector B* mesons as 
explicit degrees besides the B and light pseudoscalar <p mesons since the B and B* mesons would form a degenerate 
doublet in the limit of heavy quark symmetry. We count the mass difference A between B and B* mesons as O^e 1 ). 
The potentials of the BB mesons starts at O(e ). We will investigate the corrections up to 0(e 2 ). 

This paper is organized as follows. In Sec. [TTJ we list the Lagrangians of HMxPT. In Sec. |III| we present the 
expressions of the BB potenti al, w hich include the tree-level diagram contributions at the leading order and the loop 



corrections at 0(e 2 ). In Sec. IV we give the numerical results of the BB potentials in the first subsection. Then 
we present the results of the potentials of the DD mesons in the second subsection. We compare the results within 
different schemes in Sec. |Vj Sec. [VlJ is a short summary. 

II. LAGRANGIANS WITH HM^PT 

The leading order BB potential is at O(e ) and receives only the contribution from the tree-level diagrams made 



up of the vertices of the leading order Lagrangians in Eq. pi). The corrections start at 0(e 2 ). They contain the 
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contributions of the 1-loop diagrams generated by the leading Lagrangians and the contributions of tree diagrams 
generated by the Lagrangians at higher order. 
The leading Lagrangians are 

cf H = D a Tr[H^H]Tr[H^H] + D b Tr[H lfll5 H]Tr[H^ l5 H] 

+E a Tr[H^X a H]TT[H^X a H] + E b Tr[H W X a H]Tr[H^ l5 X a H], (3) 

- -<(*«• dH)H) + (HvTH) + g(H fa 5 H) -^A(Ha^Ha^), (4) 

where the number in the round bracket represents the chiral dimension, = (1, 0) is the velocity of a slowly moving 
heavy meson, and H represents the B and B* doublet in the heavy quark symmetry limit, 

H = ^{p;^+iP l5 ), ^ = 7 o jff t 7 o = ( P; t 7 M +i pt 75 )l+i ) (5) 

P=(B-,B Q ,B° S ), P; = (B*-,B*°,B* S ),. (6) 
The pseudoscalar meson field, chiral connection and axial vector field are defined as follows, 

^ = ^M, u»= l -{?,d^}, £ = exp(#/2/), (7) 
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The Lagrangian 4j? generates the contact interaction terms of the four bottomed mesons while depicts the 
interaction between the heavy mesons and light pseudoscalar mesons. The other contact terms with different Lorentz 
structures at the leading order are not independent. Actually they are linear combinations of terms in 4ff> so we 
do not need them. For example, the term TrlHj^HHj^H] can be expressed as the linear combination of terms in 

4# by Fierz transformation. In the heavy meson limit, Tr[H H]Tr[H H] and Tv[Ha tlv H] r £r[Ha p,u H] can be absorbed 

by readjusting the coefficients D a and D b of 4ff> respectively. The term Tr[Hj 5 H]Tr[H"f 5 H] vanishes in the heavy 
meson limit. Similar conclusions hold for the terms containing A a such as Tr[H-f5X a H]Tr[H-f 5 X a H] etc. 

The ranges of the couplings D a , D b , E a , and E b were estimated in the DD case by fixing the mass of X(3872), 
X(3915), and F(4140) and the isospin breaking branching ratio of X(3872) in Ref. [3T]. Their values lie from several 
to tens of GeV -2 with positive signs. 

The 0(e 2 ) Lagrangian d\ H will also contribute to the potentials, which reads 

4^ = D':Tv[H^H}Tt[H^H}Tv( x+ ) + D^ r Tr[H W H}Tr[H^ l5 H]Tr( X+ ) 

+E^TT[H 1 ^X a H]TT[HrKH}TT( X +) + E' b l Tv[H W X a H}TT[Hr l5 X a H}TT(x + ) (9) 
cf H d) = D d a Tt[H^x+H]Tt[H^H] + DiTr[H WX +H]Tr[H^ 75 H] 

+E *d abc Tr[H^X a H}TT{H^X b H}Tr\x + X c ] + ^d afcc Tr[i/ 7p75 A aJ ff]Tr[i/r 75 A 6j ff]Tr[x+A c ] (10) 
4h ) = {^aiTr[(« • DH) lfl (v ■ DH)]Tr[H^H] + Dl 2 Ti[(v ■ DH)^H]Tr[(v ■ DH)-fB\ 

+D v a3 Tr[(v ■ DH)^H]Tr[H^(v ■ DH)] + D v aA Tr[((v ■ D) 2 H)^H]Tt[H^H] 

+D v bl Tr[(v ■ DH)w 5 (v ■ DH)]Tr[H^ l5 H] + ... + E^Tr[(v ■ DH) ltl X a (v ■ DH)]Tv[H^ X a H] + ... 
+E^Tr[(v ■ DH)^ l5 X a (v ■ DH)]Tr[H^ 5 X a H] + ...} + H.c, (11) 

4h 9) = {D^[(D»H) W (D u H)]Tr[Hw 5 H} + D^[{D^ H) W H]Tr[{D v H) lv ^H] 
+DITy[{D»H) 1 ^H]T?[H 1v15 {D v H)] + DlTr[{D"D l 'H) W H]Tr[H 7v75 H] 

+ElTrl(D»H) %ll5 X a (D»H)]Tr[H lu "f 5 X a H] + ...} + H.c, (12) 

where d abc is the totally symmetric structure constant of SU(3) group, and 

X±=X±~ ^Tr[x±], x± = C f xC f ± £x£, X = diag(m 2 , to 2 , 2m 2 K - m 2 ). (13) 
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The low-energy constants (LECs) in Eqs. ([£])-( 12 ) contain both the infinite and finite parts. We will use the infinite 
parts to cancel the divergence introduced by the loop diagrams. We are unable to determine the finite parts of the 
low-energy constants (LECs) due to lack of experimental data right now, which we tend to neglect for the moment. 
However, these LECs at 0(e 2 ) should be included in a complete analysis in the future when more experimental data 
are available. 

^ah an d ^iff are ma de up of four heavy mesons fields, Tr(x+) and traceless \+- The LECs in C^2h^ anc ^ QiiP 
will absorb the divergent parts from the one- loop diagrams which are proportional to m 2 . will absorb the 

divergent parts proportional to the square of the external line energy. There are also divergent parts proportional to 
the mass difference A, which will vanish in the heavy meson symmetry limit. These divergence can be absorbed by the 
additional four heavy meson interaction terms proportional to A. does not contribute to the renormalization 

of the BB potentials. Instead it will contribute to the BB* and B*B* potentials. 



III. POTENTIALS WITH HMxPT 

With the strict isospin symmetry, there are only four independent potentials for the channels BB , BB°, B S B S °, 

- - 1/2 

and BB S . The superscript represents the isospin of the channel. At the leading order, the potentials of the 
bottomed mesons only receive the contributions from the tree diagrams with the contact terms in C^, 

ySl=-^-\ E ^ V S%=°> V ( £ Bs0 =-lD a -\E a , V^=-2D a -lE a . (14) 

The loop diagrams in Figs. [2] and [3] will contribute at the next to leading order. The diagrams h.l, h.2, and B.l in 
Figs. [2] and [3] contain both 2PR and 2PI parts if using the ordinary Feynman rules of HM;\PT. We must remove the 
2PR contribution to get the correct potentials. The 2PR parts result from the double poles of the two heavy mesons, 
which can be removed by the careful subtraction in the propagator of the heavy mesons 

1 



[v ■ pi + Si + is] [v ■ p 2 + 82 + is] 



— 2ni5(v ■ pi + Si)] — > — 2 v ■ P2 + S 2 = —v ■ pi — Si 



v ■ pi + Si + is v ■ pi + Si + is (v ■ P1 + S1 + is) 

1 



(15) 



[v ■ pi + Si + is] [v ■ p 2 + S 2 + is] 



other 



We calculate these Feynman diagrams with dimension regularization and modified minimal subtraction scheme. 
The divergent terms proportional to L will be absorbed by the contact terms at 0(e 2 ), where 

L= T \ + -(7 B -l-ln4ir) . (16) 

Here is the Euler constant 0.5772157, A is the scale of the dimension regularization, and we set A = Airf. 








g-1 g-2 h.l h.2 z.l 

FIG. 2: The loop diagrams with a contact vertex. The thin solid, thick solid and dashed lines represent the heavy pseudoscalar 
mesons, heavy vector mesons, and light pseudoscalar mesons respectively. 

The potentials are finite after the renormalization of the wavefunctions and vertices. The diagram z. 1 in Fig. 
arises from the renormalization of the wavefuntions. The combined divergence generated by diagrams g.l and g.2 in 
Fig. [2] can be absorbed by the LECs D^, E v ai , D^ d ', and E^ d . The divergence generated by the diagram h.l or h.2 
can be absorbed by the redefinitions of LECs in cf^\ C±ff\ and cf^\ For the 20 exchange diagrams in Fig. [3J 
the divergence of the football or triangle diagram can be absorbed by E^, and E^ d '. The divergence of the box and 
crossed diagram can be absorbed by D^, D^ d , E v ai , and E^ d . 
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F.l T.l T.2 B.l R.l 



FIG. 3: The 2</>-exchange diagrams including the football diagram (F.l), triangle diagrams (T.l and T.2), box diagram (B.l), 
and crossed diagram (R.l). 



The potentials obtained from the same Feynman diagram for the different channels differ just by a flavor dependent 
coefficient. So it is convenient to write down the potential of the channel l ch' in the form 



r 



ch 



E 

diag.mi ,m,2 , 



. h (m 1 ,m 2 , ...) ® Y diag {m u m 2 , ...), 



(17) 



where '<iia<?' runs over all the Feynman diagrams shown in Figs, ^andpl and runs over {m w ,m^-,m^}. Y is a 
scalar function independent of the flavor structure of the channel 'cV, while (3 is the flavor dependent coefficient. 
The corresponding Y functions of Fig. [2] are 



Y gl (m 
Y g - 2 (m 
Y hl (m 
Y h - 2 (m 
Y zA (m 
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•T ~ 

D 1 
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]J| 2 | (m,£-«,-A,£-A), 

]J| 2 | (m,£ + ®,-A,£-A), 

]j£ 2 } (m,£-A,£-A), 

hj^] (m,£-q -A,£ + q -A), 

r 

]d x J$ 2 \ {m,x) 



(18) 
(19) 
(20) 
(21) 
(22) 



where we work in the center-of-mass frame of the incoming heavy mesons, £ is the residual energy of the incoming 
heavy meson (the difference between the energy and the B meson mass), and q is the transferred momentum. The 
definitions of the J functions are collected in Appendix [Aj {X} r represents the finite part of X , 



\X} r = lim (X - L-S-X) + lim (JLJLx). 



(23) 



The loop diagrams in Fig. [2] are made up of the contact vertices of £4^, so the flavor dependent coefficient /3 can 
be written down as 



= DJ 



(24) 



Actually, only three coefficients /3s for the diagrams g.l, h.l, and z.l in Fig[2]are independent. For convenience, we 
list {^D^^Da^Db^Db} in Tableland {/3£ a \ /3£ a \ /3| b \ ^} in Table [n] The coefficients /3s for g.2 and h.2 can be 
obtained by the following relations 



(25) 



59-2 _ og.l og.2 _ og.l oh.2 _ 5ft. 1 og.2 _ og.l og.2 _ og.l oh. 2 _ 2ft. 1 
PDa — PDa > PDb — PDb ' PDb ~ PDb ) PEa — PEa ' PEb — PEb ' PEb ~ PEb > 

and all the others are zero. 

Now one can write down the potentials induced by Fig. [2] for the different channels. Let us take the BB° system 
as an example. One can get the contribution from the diagram g.l by Tables [T| and |TT| 



VK-L = - l -[2AD a -2AD b + i2E a -'62E b ]Y3 1 (m 7T )- 1 - 



BB" 



32 

16D a ~Y E -~ 32Eb 



TABLE I: The coefficients /3s of the loop diagrams with a contact term: {/3^ a , j3^ a , /3^ b , /3 Db } 
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TABLE II: The coefficients /3s of the loop diagrams with a contact term: {/3^, (3^, Pkb} 
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The diagrams in Fig. [^represent the potentials with the 2<j> exchange. We list the /3 functions {/3 s ' 1 , ft RA } in 
ITi*} and all the others can be obtained by the following relations 



The corresponding Y functions of Fig. [3] are 
1 
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+ hi^JS + h-Df ~?]J& + [-~D? V] J« + [-Af?]J« 



+[\{??]J§i + [\(?)V& \ (">■ v. e a. e - a. 7) . 



Again, taking the BB° channel as an example, the potential from the diagram B.l of Fig. [3]reads 



V BB° ~ 



-QY^m*,™,) - -Y* L (m n ,m ri ) + Y* A (m ni m n ) +Y aA {m ni m 7 <) 
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TABLE III: The coefficients /3s of the 2-0 exchange diagrams: {/3 S1 , 



(m w , rn„) (m K , m K ) K, m^) (m„, m K ) (rn^, m n ) (m K , m n ) (m K , m^) (m,,, m„) (m^, itik) 
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Finally, the potentials at 0(e 2 ) can be obtained by summing the products of the corresponding /3 and Y as in 
Eq. Jl7b. ' 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



Potentials of BB Mesons 



We have calculated the potentials of BB mesons up to 0(e 2 ) for four independent channels. The 0(e 2 ) potentials 
contain two parts 

y-(2,cont) and y(2,20)^ corre sponding to Fig. § and Fig. [3] respectively. We will focus on the 
potentials with £ = q° = in the momentum space. After the Fourier transformation, we can get the traditional 
potentials in the coordinate space. The other parameters are listed as follows [32H34] 



139 MeV, m K = 494 MeV, 



ml) A = 46 MeV, / = 92 MeV, g = 0.52, (34) 



where the A is the mass difference between B and B* mesons, / is the pion decay constant, and g is the axial coupling 
constant from the unquenched lattice QCD simulation. 

The potentials and T/( 2 > cont ) are both generated by the contact vertices. They are independent of the transferred 
momentum \q\. They are (5(r)-like potentials in the coordinate space. From Eq. (14), we notice that the terms 

proportional to Db and Eb in do not contribute to V^°\ The potential vanishes in the channel (BB) at the 
leading order. At the next to leading order, the situation is different: 



V, 



(2,cont) 



BB 1 



-0.32.E Q - 0.32E b , V. 



(2,cont) 



BB 



0.19A, - 0.02L> b - 0.085£ a - 0.36E b , 



V BB™ t] = -0-53S a - 0.53£ b , 



V BS$ = -°- 32 ^a - 0.32£; 6 . 



(35) 



All the terms in contribute to y( 2 > cont ) . However, in the BB 



1/2 

B S B S , and BB S channels, the contributions 
proportional to D a or Db from different diagrams in Fig. [2] cancel each other. Roughly speaking, the corrections 
are small compared with the leading order contribution. We also notice that 



(2,cont) 



B S B S 



> 



V 



(2,cont) 



BB 



1/2 



V, 



(2,cont) 



BB 



(36) 



As we have emphasized in the previous section, the finite parts of the 0(e 2 ) LECs in Eq. ([9|-(12| also contribute 
to the potential while its divergent parts cancel the divergence from the one-loop diagrams. Unfortunately we are 
unable to fix these LECs because of lack of experimental data. In the following analysis, we focus on the behavior of 
the 20- exchange potentials. 

We plot the 20-exchange potentials V^ 2,2 ^ of the BB mesons in Fig. [4} From the figure, the contributions from 
the football and triangle diagrams are coincidentally close in all the BB channels. When the transferred momentum 
is small, the sign of the potential from the crossed diagrams is different from the other 2</>-exchange diagrams'. We 

have noticed that the 2</>-exchange potentials of the BB 1 , B S B S °, and BB^ 2 channels are negative in the small- 
momentum region. In other words, the 20-exchange interaction is attractive if we ignore contribution from the LECs. 
In contrast, the 20-exchange interaction of the BB° channel is repulsive in the small-momentum region without these 
LECs. The 20-exchange potential in the B S B S ° channel is nearly twice larger than those in the other channels. 

The potentials from the football diagr am, t rian gle diagram, box diagram, and crossed diagram are proportional to 
3°, <7 2 , <? 4 , an d g A respectively from Eqs. ( 28 )-( 32 ) . So due to the coupling constant g = 0.52 is quite small, one would 
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FIG. 4: The BB potentials with the 2<f> exchange. 



naively expect that the potential from the triangle diagram is suppressed by a factor about 0.27, and the potential 
from the box or cross diagram is suppressed by about 0.07 compared with that from the football diagram. However, 
we do not see the suppression in the Fig. |4| There is also an enhancement of the flavor coefficient f3 for the triangle, 
box, and crossed diagram from Eq. ( |27[ ), which roughly compensates the suppression of the small g. That is why we 
sec neither the suppression due to the small g nor the enhancement of j3 in the numerical results. If we let g — > 1, the 
potential from the box or crossed diagram would be much larger than the potential from the triangle diagram, which 
would be larger than the potential from the football diagram. 

From Fig. |ZJ we notice that the contribution from the box diagram dominates the potential V 1 - 2 ' 2 ^ . From Eq. 



(32 1, we have 



15 < y R1 K,^ 3 ) < 50; (37) 

~ Y BS (m w ,in n ) ~ 

where the intermediate meson pair (frifa can be ttK, ttj], KK, Krj 1 or rjrj. So the potential from the box diagram is 
dominated by the intermediate states with at least one kaon or eta meson. 



B. Potentials of the DD Mesons 



Similarly we can study the potentials between the D°, D + , and mesons. Now the intermediate heavy vector 
mesons are D*°, D* + , and D* + . The mass difference A increases to 142 MeV. The axial coupling constant g = 0.59 
from the decay width of D* + [32]. The LECs D a , E a etc should be modified correspondingly. The expressions for the 
DD mesons are the same as those for the B5's except that the channels are DD 1 , DD°, D S D S °, and DD 1 ! 2 . Thc 
DD potentials with the 20 exchange are plotted in Fig. [5] The potentials related to the contact terms are 

^2,cont) = QmE DD + Q 89E DD^ ^co„t) = _ QMD DD + 2D DD _ 0M5 Q E DD + U£ DD 

Vg£f = 1 . IE° D + 1 .IE° D , V%£*$ = 1.1E° D + 1 .1E° D . (38) 
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FIG. 5: The DD potentials with the 2<f> exchange. 



We notice that there is big difference between the potentials V^ 2,0011 *-* of the BB and DD mesons by comparing 
Eqs. (35) and (38). The difference originates from the different axial coupling g and mass difference A. The signs 



of most of the terms in V^ 2 ' cont ^ are different for the bottom and charm cases if assuming D^ (E^f) is equal to 
D a /b(E a /b) for the B mesons. One obtains the relation 



V, 



(2,cont) 



V, 



(2,cont) 



D S D S 



> 



V, 



(2,cont) 



DD 1 



(39) 



which is different from that in the case of the B mesons. 

Comparing Fig. 4 and Fig. pi one can find that the total Vgg^ is roughly twice of V D 2 D 2 ^ for each channel. 
Moreover, the separate contributions from the crossed or triangle diagrams have opposite signs for the B and D 
mesons. 

One can obtain the potentials of anti-heavy mesons based on C-parity conservation 



V, 



BB 1 



BB 1 



v, 



DD 1 



DD 1 



(40) 



V. COMPARISON BETWEEN RESULTS IN DIFFERENT SCHEMES 



One can also systematically study the potentials of the heavy pseudoscalar mesons without heavy vector mesons as 
the explicit degrees. The contributions from the heavy vector mesons, as well as other resonances, will be embodied 
in the low-energy constants. In the scheme without heavy vector mesons, the potential at the leading order remains 
the same. However, only the football diagram survives at 0(e 2 )- 

It is also very interesting to investigate the case with strict heavy quark spin symmetry. Now the heavy vector 
mesons are included as the explicit degrees but the mass difference A is set to be zero. When A approaches 0, the 
potentials induced by the diagrams h.l, h.2 and B.l will approach infinity if the 2PR contributions are not removed. 
The two-heavy-vector-meson-reducible contribution does not appear when the mass difference A is finite. It only 
appears as A = before the subtraction. So there is a jump in the potentials as A goes from nonzero to zero. When 
solving the nonperturbative equation to get the observable, the potentials such as Vbb^B'B* m ust be included as 
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A = 0. But they are not necessary as A ^ 0. Thus the jump might vanish for the observable such as the binding 
energy. 

In the new approximations, the potentials do not change at the leading order. The difference appears at the next 
to leading order. In the approximation without the heavy vector mesons, there does not exist those diagrams in Fig. 
[2] so y( 2 - cont ) = 0. In the approximation including the heavy vector mesons and A = 0, we have 

y(2,cont) = _ .78£: Q - 0.78£ 6 , Vl 2 - C ont) = 0.48L> Q - 0.077A, - 0A7E a - 0.91E b , 

V tj}? = ~ h2Ea ~ L2 ^ fc ' V SB™> = -0-8SE a - 0.88£ b , (41) 

y [ DD° nt) = ~ l - QE a ° - 1-0E? D , Vg£° nt) = 0-62D° D - 0.099£ b DD - 0.22E? D - 1.2E? D , 
V d1d? = ~ 1 - QE a° - l - &E b° > Vg£°$ = - 1 . LE? D - 1 . IE? . (42) 
The difference between the BB and DD potentials only originates from the different axial coupling g as A = 0, 

F (2,cont) 2 

_J3B = = 0.8, forA = 0. (43) 

Tr(2,cont) Q.59 2 V ' 

V DD 

The BB potential with A = is nearly twice as large as that with A = 46 MeV in every channel by comparing 



and 



Eqs. (35) and (41). The potentials with A — 46 MeV are approximately equal in the channels BB and BB S 



while those with A = are not equal. The difference between the DD potentials is even larger with different A from 



Eqs. (38) and (42). The sign of the DD potential as A = is different from that as A = 142 MeV in every channel. 
The potential of the channel D S D S ° is very close to that of the channel DD^ 2 for the case with A = 142 MeV, but 
the situation is different for the case with A = 0. 

The difference between the potentials with different A mainly results from the subtraction of the 2-heavy-vector- 
meson-reducible contributions to get the potentials of the heavy pseudoscalar mesons as A = 0. To recover the 
2-heavy-vector-meson-reducible contributions, one should include the potentials such as Vgg^g.g, when solving the 
nonperturbative equations to get the observable as A = 0. 

We list the BB and DD potentials with 2<p exchange in different approximations in Figs. [6] and [7] respectively. 
From Fig. [6j the potentials are relatively close between the case without B* (case I) and that with A = (case II). 
The potential |V| for case I or II is about 15%-50% of that for the case with B* and A = 46 MeV (case III) in every 
channel. The dependence of the potential on \q\ for case III is slightly stronger than that for the other two cases in 

the channels BB 1 , B S B S ° , and BBg 1 ' . The potential decreases for case III and increases for the other two cases as 
\q\ grows in the channel BB° . 

The situation in Fig. [7] is similar to that in Fig. [6] But the difference of the DD potentials between different cases 
is smaller than that of the BB potentials. We notice that the DD potential is equal to the BB potential for case I 
in each channel since the potential from the football diagram is independent of the mass difference A and the axial 
coupling constant g. 

VI. SUMMARY 

In a short summary, we have calculated the potentials of the heavy pseudoscalar mesons up to 0(e 2 ) in the 
momentum space with HM^PT. We have carefully analyzed the tree-level contribution and one-loop correction to the 
contact vertices, and the 27r-exchange contribution. We have also discussed the potentials in different schemes. 

Generally speaking, the potential of hadrons can be separated into the long-range, medium-range, and short-range 
parts. For the two heavy pseudoscalar mesons, there does not exist the long-range l(/)-exchange potential. The 
medium-range potential contains the 2^-exchange potential and the contributions by the Lagrangians ®-{l2). The 



20-exchange potential is model-independent since there are no unknown constants in it, which is very essential for 
the medium-range interaction of the heavy pseudoscalar mesons. The interaction induced only by the 2<p exchange 
is repulsive in the channels BB 1-0 , DD 1-0 , while attractive in the other channels. Unfortunately the leading order 
coupling constants from the contact terms and LECs at 0(e 2 ) remain undetermined due to lack of experimental data. 

Once these LECs are extracted from lattice QCD simulation, other model approaches or future experimental 
measurements, the potentials derived in this work can be used to study the possible molecular states or scattering 
phase shift of the two heavy pseudoscalar mesons system. On the other hand, the analytical chiral structures of the 
potentials of the heavy meson pair may be useful in the extrapolation of the heavy meson interaction from lattice 
QCD simulation. 
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FIG. 6: The BB potentials with the 2<p exchange in different schemes. 
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Appendix A: Some functions used for potentials 

The J functions can be obtained by calculating the following integrals in D dimensions 



(2tt) D (v -1+uj + ie) (I 2 - m 2 + is) 
= {J a , v a J&, v^J^+g^J^ (g\/v)J^+v a v^v''J^}(m,co), (Al) 



f d u lX A ~ L 
7 72^ 



(2tt) D (I 2 -m 2 + ie) 

{J C , 0, g a0 J^}(m), (A2) 
d D ix i - D {l, i a , l a lP, i a iH~<} 



(2ir) D (v -1+UJ + ie)[(+/-)v ■ 1 + 5 + ie] s {l 2 - m 2 + ie) 
{J°'\ v a Jf( h , v^Jl( h +g^J!t, ( 9 Vv)J§l h + v a vf'v''JZl h }(m,u,,S), (A3) 

D lx i-D | 1; ^ l a lPp} 



. r d u ix 4 

J (2n) 



(2ir) D (I 2 — m 2 + ie)[(q + I) 2 — M 2 + ie] 
{j<f, fj*, q a q r j J 2 F 1+ g^j[ 2 , (gVq)j( 1+ q a q P q ^j( 2 }(m,M, q ), (A4) 

d D ix i - D {i, i a , i a iP, i a iPp, ri^pi 5 } 



(2ir) D (v-1+lu + ie)(l 2 - m 2 + ie)[(q + I) 2 - M 2 + ie] 
{J T , q a J? 1+ v a J? 2 , g^J 2 T 1 +q a q f) J 2 T 2 + v a v^J 2 T 3 + (q\/v)J 2 T i , 
(g V q )jT + q a q S\Jl 2 + ( g 2 V v) J 3 T 3 + (g V v) J 3 T 4 + ( 9 V « 2 ) J 3 T 5 + v a v^ J 3 T 6 , 

(5 V g) J 4 T ! + ( 5 V g 2 ) J 4 T 2 + q a q q 1 q s Jj i + (.g V « 2 ) J 4 T 4 + v a v^v^v 5 J 4 T 5 + (<? 3 V w) J 4 T 6 + (g 2 V v 2 )jj 7 + (gV « 3 ) J 4 T 8 
+( 9 V g VD)J 4 T 9 }(m,M, W , g ), (A5) 

d D ix 4 - D {i, p, z q P, ppt, i a i?pi s } 



(2n) D (v-1+uj + ie)[(+/-)v -1 + 6 + ie] s (l 2 - m 2 + ie)[(q + I) 2 - M 2 + ie] 
= {4 /B , q a J* /B +v a J* /B , g^J^ B + q a q f3 J% B + v a v?J 2 R 3 /B + (q V v) J?/ B , 

(.9 V q) J% B + q a q fi q 1 J*l B + {<? V v) J* /B + (gV v)J* /B + (gV v 2 ) J*' B + v a v^ jf' B , 

(9 V g)jT + (g V q 2 ) JT + tfq?q'<fj%( B + (g V v 2 ) J% B + „<»,, Wj* /B + ( g 3 v V )J$ B + (q 2 V v 2 )jf' B 
+(q V v 3 )J* /B + (gVqV v)jg /B } (m, M, u, 6, q), (A6) 

where we have used the following Feynman rule for propagators of two heavy mesons in HMxPT to remove the 2PR 
contributions from the Feynman diagrams 

1 



(v ■ I + oj + ie)[(sgn)v ■ I + S + ie] s 



' "[(-) — , 1 , • -2m5(vl + u)] B = (-) 7 . , 1 , . , 2 sgn: 



(v-l + to + ie) v-l + us + ie ' ' (v ■ I + uj + ie) 2 {A.7) 

other 



(v ■ I + uj + ie) [(sgn)v ■ I + S + ie] 
The notation X V Y V Z V ... represents the symmetrized tensor of X a YP Z 1 and in detail, 
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q 2 V V = q f3 q 1 v a + q a q l v fi + q a q P v~«, q V v 2 = q<v a v fi + q f3 v a v"< + q a v f3 v"< , 

5 V g = g al3 g' fS + g aS g^ + g al g PS , g V q 2 = q a q f3 g lS + q a q S g M + q a q 1 g fiS + q' f q S g aP + AV 7 + A 7 *?" 5 , 
gVv 2 = v a v f3 g~ y5 + v a v S g^ + v a v~<g f3S + v~<v s g al3 + v p v s g ari + v^v^g" 5 , 

q 3 Vv = qPq~<q S v a + q a q^q S v P + q a q !3 q 5 v'< + q a q f3 q'<v 8 , qW 3 = gV^V + q~ r v a v f3 v 5 + q f3 v a v~<v S + q a v f3 v~<v\ 
q 2 Vv 2 = q*q 5 v a v 13 + q f3 q S v a v' y + q a q S v l3 v' 1 + q' 3 q 1 v a v S + q a q*v' 3 v & + q a q' 3 v^v S , 

gVqVv = qPv a g~< 5 + q a v l3 g'< S + q 5 v a g^ + q<v a g l3S + q a v S g f3 ~< + q a v'<g' 3S + q S v~>g af! + q S v !3 g a '< + q'<v S g a ^ 

+q''v l3 g aS + AV 7 + q' S v^g aS . (A8) 
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